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effects  of  a  discretization  carried  out  in  the  partial  differential 
equation  for  different  formulations  of  the  far  field  conditions. 
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but  the  relations  obtained  for  different  formulations  arise  from 
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approach  an  approximation  of  the  potential  and  its  normal  deriva¬ 
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o their— words'  the  vector  is  restricted  to  finite  dimensional 
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field  conditions  depend  upon  the  choice  of  the  subspace  and  of  the 
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For  vectors  in  an  infinite  vector  space  they  are  independent  of 
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representation  for  the  potential. 
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SECTION  1 
INTRODUCTION 

In  Reference  1,  far  field  conditions  for  small  periodic 
oscillations  superimposed  on  a  subsonic  flow  have  been  derived. 

They  express  the  effect  of  the  distant  field  on  the  portion  of  the 
flow  field  that  is  to  be  treated  numerically.  One  obtains  linear 
relations  in  which  the  potential  and  its  normal  derivative  at  all 
points  of  the  far  boundary  of  the  computed  part  of  the  flow  field 
appear  simultaneously.  We  speak  of  global  conditions  since  all 
boundary  points  are  involved.  Approximate  local  far  field 
conditions  have  been  derived  by  Bayliss,  Gunzberger  and  Turker 
(Reference  2) . 

The  present  report  clarifies  by  means  of  an  example  some  of 
the  questions  which  arise  in  the  application  of  these  conditions; 
of  particular  interest  are  the  errors  introduced  by  the  discretiza¬ 
tions  which  are  necessary  in  a  numerical  treatment.  The  example 
is  rather  simple;  it  can  be  treated  analytically.  The  analytical 
properties  will  be  helpful  in  the  discussion  but,  since  one  wants 
to  imitate  more  realistic  problems,  they  will  not  be  used  in  a 
numerical  approach.  Actually,  the  numerical  solutions  play  a  very 
minor  role.  The  problem  has  been  programmed  and  solved  for  different 
boundary  conditions ,  but  only  to  check  the  feasibility  of  the  approach, 
not  to  provide  a  production  routine. 

A  comparison  of  the  results  so  obtained  with  exact  solutions 
will  not  be  carried  out  because  it  does  not  allow  us  to  separate 
different  error  sources  (errors  due  to  the  discretization  in  the 
differential  equations,  in  the  far  field  conditions,  and  in 
satisfactory  boundary  conditions  at  the  profile).  Instead,  we 
shall  examine  how  the  discretization  affects  different  eigenfunctions. 

In  realistic  problems  the  differential  equation  governing 
unsteady  perturbations  is  not  identical  with  the  one  for  perturbations 
in  a  parallel  subsonic  flow,  in  particular  not  in  the  vicinity  of  the 
profile.  This  does  not  limit  the  scope  of  the  present  study  because 
in  the  distant  field,  the  latter  gives  an  acceptable  approximation. 
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SECTION  II 

SOME  PRELIMINARY  OBSERVATIONS 


We  consider  the  unsteady  two  dimensional  field  generated  by 
a  plate  which  extends  along  the  x  axis  of  a  Cartesian  system  of 
coordinates  x,y  from  x  =  -1  to  x  =  +1.  Along  the  plate  the  velo¬ 
city  normal  to  it,  that  is  the  normal  derivative  of  the  oscillating 
part  of  the  potential  <)>  is  prescribed.  For  an  infinitely  thin 
oscillating  plane  the  value  of  <j>  at  its  upper  and  lower  side  are, 
of  course,  the  same,  but  we  shall  allow  for  a  more  general  situation. 
The  problem  has  a  number  of  symmetries  which  allow  one  to  separate 
it  into  smaller  problems;  in  the  subproblem  treated  numerically  here 
the  potential  <J>  is  symmetric  with  respect  to  the  x  axis.  The 
components  of  the  perturbation  velocity  then  have  opposite  signs  at 
the  upper  and  lower  side  of  the  plate.  To  comolete  the  formulation 
of  the  boundary  value  problem  one  must  prescribe  at  a  distance  the 
"far  field"  conditions  under  study  in  the  present  report. 

The  equations  for  the  distant  field  can  ultimately  be  reduced 
to  the  Helmholtz  equation 


b  +  <b  +  p  <}>  =  o 
xx  Yyy  v 


(1) 


It  arises  from  the  eauation  of  two  dimensional  acoustics 


<j>  + 

xx 


b  -  (J\  .  =  0 
YY  tt 


(2) 


by  setting 


<j)  (x ,  y ,  t )  =  <{>(x,y)  exp(iyt)  (3) 

In  this  formulation  it  is  assumed  that  the  time  has  been  scaled 
so  that  the  velocity  of  sound  is  1.  The  time  for  one  period  of 
the  oscillatory  motion  is  then  T  =  2ir/u ;  this  formula  gives  also 
the  wave  length  (because  of  the  choice  of  the  velocity  of  sound) . 


To  get  some  appreciation  for  the  magnitude  of  it  we  mention  that 
its  value  is  tt  for  a  wave  with  a  length  equal  to  the  width  of  the 
plate  (which  is  2  in  the  present  example) . 


If  the  plate  is  embedded  in  a  subsonic  parallel  flow,  one 
must  make  several  transformations  in  order  to  arrive  at  the 
Helmholtz  equation  (see  Ref.  1).  If  in  such  a  flow  the  plate 
oscillates  with  a  circular  frequency,  v,  then  one  obtains  (after 
the  transformation  to  the  Helmholtz  equation) 


aL  _  vL 

2  2  _  2 
-U  a(l-M 


where  U  is  the  free  stream  velocity,  "a"  is  the  free  stream 
velocity  of  sound,  M  =  U/a  is  the  free  stream  Mach  number,  and  L 
is  some  characteristic  length  (here  the  half  chord  of  the  plate) . 

During  one  period  T  =  27t/v  a  perturbation  travels  strictly 
upstream  or  strictly  downstream  by  the  respective  wave  lengths, 
via  (  2tt/v)  a  (1-M)  and  ( 2u/v)  a  (1+M)  . 

Frequently  a  reduced  frequency,  here  denoted  by  u 1 ,  given 
by 


y '  =  v  L/  a 


(5) 


is  introduced.  Then 


M 

y  =  y  — — T  (6) 

1-M 

The  frequencies  v  are  determined  by  the  physical  situation 
under  investigation. 

The  transf ormation  to  the  Helmholtz  equation  generates  a 
wave  pattern  in  which  the  wave  length  is  the  same  in  all  directions. 
For  a  given  v  one  has  the  following  expressions  for  the  wave  length 


Helmholtz  equation 


2v,/a  =  ( 2 r/vL) a ( 1-M2 ) 


waves  travelina  upstream 
in  a  parallel  flow 


(?■ 


1.)  a  ( 1-M)  =  ( 2  "/\  L)  a  ( 1 


■M2) (1+M) _1 


wav^o  ui'aveling  downstream 
in  a  parallel  flow 


(2i:/vl)  a  (1+M)=  (2  -i/vL)  a (1-M  )  (1-M) 


The  wave  length  for  the  Helmholtz  equation  is  the  harmonic  mean 
between  the  wave  lengths  for  upstream  and  downstream  traveling 
waves . 

The  fact  that  after  the  transformation  to  the  Helmholtz 
equation  the  wave  length  is  larger  than  for  waves  traveling  upstream 
in  the  original  field  is  of  practical  interest.  The  grid  points 
of  a  finite  difference  net  must  be  close  enough  to  reproduce  the 
waviness  of  the  unsteady  field.  After  the  transformation  to  the 
Helmholtz  equation  a  grid  which  is  coarser  than  that  of  the 
original  problem  is  admissible.  Of  course,  for  small  values  of  v 
and  Mach  numbers  not  too  close  to  one  this  is  less  important,  for 
the  wave  lengths  are  large  to  begin  with  and  the  solution  resembles 
in  essence  that  of  the  Laplace  equation.  But,  under  different 
circumstances,  it  may  be  advantageous  to  carry  out  such  a 
transformation  even  though  the  differential  equation  for  the  vicinity 
of  the  profile  differs  from  the  Helmholtz  equation. 
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SECTION  III 
TRANSFORMATIONS 


In  a  numerical  approach  one  must  provide  a  net  in  which  the 
discretization  is  carried  out.  One  might  use  a  rectangular  net 
in  the  physical  plane.  An  elimination  strategy  for  this  purpose 
is  described  in  Reference  3.  In  Reference  1  a  different  net 
obtained  by  a  conformal  mapping  has  been  suggested.  For  practical 
purposes  the  first  net  is  probably  preferable  because  it  retains 
the  orientation  imposed  by  the  free  stream  direction.  The  net 
of  Reference  1,  used  here,  makes  it  possible  to  study  the  problem 
in  greater  detail.  The  general  conclusions  do  not  depend  upon  the 
choice  of  the  net. 

The  starting  point  is  the  Helmholtz  equation,  Eq .  (1). 

Let 


X  + 

iy  =  z 

(6) 

K  + 

in  =  t 

(7) 

and 


let 


z  =  g  (  ', ) 


g  '  =  dg/dj;  . 

where  g  is  an  analytic  function.  Then  one  obtains 

2  i  2 

4>r  +<j)  +]J  j  g  ’  <J>  =  0 

££  r\H  1  y 


We  choose  specifically  the  familiar  mapping 

c;  (  )  -  cosh  r. 


Then 


x  =  Re  g  =  cosh  f  cos  n 
y  =  im  g  =  sinhf  sin  n 


(8) 


(9) 


(10) 


(ID 
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for  £  =  0  one  obtains  y  =  0  and  x  varies  from  -1  to  +1  and  back 
as  n  varies  from  0  to  2 tt .  This  is  the  map  of  the  plate  surface. 
Both,  the  lines  n  =  0  and  n  =  2 tt  map  into  the  x  axis  from  x  =  1 
to  x  =  «,  the  line  n  =  n  maps  into  the  x  axis  from  x  =  -1  to 
The  whole  flow  field  is  mapped  into  a  strip  of  width  2tt  which 
extends  from  £  =  0  to  £  =  °°. 

One  obtains,  by  eliminating  n 

2  2 

_  =  ! 

2  2  x 
cosh  £  sinh  £ 

and  by  eliminating  £ 


cos  n  sin  n 


The  lines  £  =  const  and  n  =  const  are,  respectively,  ellipses  and 
hyperbolae  with  foci  at  x  =  +1  and  x  =  -1. 

A  number  of  formulae  are  listed  for  future  use.  One  has 


g'  =  sinh  (£)  (12) 

,,2  * 

| g ' I  =  (sinh  £)(sinh  L) 

where  the  star  denotes  the  conjugate  complex.  Writing  the 
hyperbolic  functions  in  terms  of  exponential  functions  and 
carrying  out  the  multiplication  one  finds 


6 


and  ultimately, 


r  =  [1/2  (cosh(2fJ  +  cos(2n)] 


1/2 


For  F,  large,  one  has 


and  therefore 


one  observes  that 


cosh  2C  ^  1/2  exp  ( 2 £) 


r  %  1/2  exp  C 


exp  ( i6 )  =  =  alii  =  ,5^) 

p  r  r  |  g  U ) 


Hence,  asymptotically  for  large  values  of  f 


exp(i0)  %  exp (in) 


We  note  that,  because  of  Eq .  (16) 

cos  (n0 )  =  Re ( exp ( i0 ) ) n  =  Re [ (g ( £ ) /r ) n] 

sin (n0 )  =  Im(exp(i0))n  =  Im [ (g U ) /r ) n] 


(16) 


(16) 


(17) 


(18) 


These  expressions  are  readily  evaluated  by  complex  arithmetic. 
Considering  x  and  y  as  functions  of  £  and  n?  one  obtains 


9x/3£  =  Re  g ' 
3x/9n  =  -im  g' 
3y/3£  =  Im  g ' 
3y/3n  =  Re  g ' 


(19) 
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Furthermore,  if  f,  and  n  are  considered  as  functions  of  x  and  y, 


3 5/ Ox  = 

|g  '  1  ~2  Re  g  ' 

oryov  = 

I  g  '  I-2  Im  g' 

3n/3x  = 

- | g ' |  ^  Im  g 

3n/3y  =  | g ' | ~2  Re  g ' 


(20) 


Since 


log  g  =  log  r  +  i0 


one  obtains  in  analogy  to  Eq.  (19) 

3  log  r/35  =  Re(g'/g) 

3  log  r/3n  =  -lm(g'/g) 
36/35  =  Im(g'/g) 

30/3n  =  Re(g'/g) 


By  substituting  Eq.  (14)  into 

g 1  _  Re  g'  +  i  Im  g* 
g  Re  g  +  i  Im  g 

one  finds 

g'/g  =  (1/r ) [ (Re (g ' )  cos  8  +  (Im(g'))sin  0 

+  i ( (Im(g ' ) cos  0  -  (Re(g')  sin  9)] 
One  has,  specifically 


g 1  sinh  ~ 
g  cosh  ’ 


(21) 


(22) 
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Introducing 


exp  C  =  u  (23) 

and  separating  the  real  and  imaginary  parts  one  obtains 


2  -2  „ . 

u  -  u  +  2i  sm  2n 

~2  ^2 

u  +  u  +2  cos  2r\ 


(24) 


Hence,  asymptotically,  for  large  u 

Re(g'/g)  =  1 

-2 

Im(g/,g)  =  2u  sin  2n 

Substituting  Eq.  (13)  into  Eq.  (9) ,  one  obtains 

<t>r^  +  +  (W2/2)  (cosh  2£  -  cos  2n)(j)  =  0 


(25) 


(26) 


Replacing  £  by  u,  (Eq.  (23)),  one  obtains  a  form  that  is  more 
convenient  for  asymptotic  considerations  and  in  addition  gives  a 
suitable  mesh  for  numerical  work 


uu 


-1, 
u  <f> 


u 


+  u 


-2 


nn 


U2l  (1+u'4)/4  - 


(u  2/2)  cos  ( 2 n >  1  4>  =  0 


:27) 


SECTION  IV 
PARTICULAR  SOLUTIONS 


The  far  field  conditions  (for  the  Helmholtz  equation),  Eq .  (57) 

of  Ref.  1,  specialized  to  the  case  of  plane  flows  is  given  by 

/lu>(x,y)<f>  (x,y)  -  a)  (x,y)  <J>  (x,y)  ]ds  =  0  (28) 

C 

Here  C  denotes  the  far  boundary  of  the  computed  part  of  the  flow 
field.  The  subscript  n  denotes  the  derivative  normal  to  this 
boundary,  ds  is  the  line  element,  <J> ,  defined  in  Eq .  (3),  is  the 
potential  expressing  harmonic  perturbations.  In  the  distant 
field  <J>  satisfies  Eq.  (1)  ,  in  the  near  field,  it  will  in  general 
satisfy  a  more  complicated  differential  equation.  The  functions  w 
have  the  character  of  test  functions,  they  satisfy  Eq.  (1)  at 
least  in  the  distant  and  usually  also  in  most  of  the  near  field. 

In  addition,  we  define  the  test  functions  used  in  the  present  report. 
They  must  satisfy  the  far  field  conditions.  One  set  is  already 
given  in  Reference  1. 


Hm2) (ur) 
m 


icos  mG 
( sin  mO 


(29) 


The  second  set  is  obtained  by  a  product  hypothesis  in  the  ' , n  plane, 
namely 


w  =  (f, ,n)  =  fk(S)gk(n) 


Introducing  a  separation  constant  one  then  obtains 


V-4  L.  , 

— i +  {~\ +  V cosh  {2r)  >  fk  =  ° 

d  r 

,  2  „ 

d  g  2 

- \  +  (:‘k  ~  y  cos  (2n))gk  =  0 

dn 


(30) 


(31) 


(32) 
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t  The  second  of  these  equations  (with  a  different  normalization  of 

|  the  constant  y2)  is  the  equation  for  Mathieu  functions,  the  first 

one  is  that  for  radial  Mathieu  functions.  One  might  refer  to  the 
literature  for  further  developments,  but  for  the  present  purposes 
it  is  nearly  as  convenient  to  derive  the  necessary  equations 
■  directly  and  to  generate  the  information  needed  by  the  computer. 

i  • 

The  expressions  (29)  and  (30)  can,  of  course,  also  be  used 
I  to  represent  the  functions  <t>  in  the  distant  field. 

|  The  periodicity  conditions  for  4>  carry  over  to  the  functions 

| 

■-  gk-  One  must,  therefore,  determine  the  separation  constants  in 

:  such  a  manner  that 


gk(°)  =  9k(2TT) 

and 

dgk/dn |  =  dgk/dn |  (33) 

n=o  n=2iT 

The  operator  of  Eq.  (32),  together  with  the  periodicity 
conditions  (33),  is  self  adjoint.  The  eigenvalues  >k  are  real,  the 
eigenfunctions  are  real  and  orthogonal  to  each  other.  The  symmetries 
of  the  problems  suggest  the  following  hypotheses 


00 


The  individual  functions  occurring  in  the  right-hand  sides  of 

Eq.  (34),  including  the  constant  term  in  the  first  expression, 

have  the  same  normalization  constant  it  over  the  interval  0  <  n  <  2it. 

The  expressions  Eqs.  (34)  are  substituted  into  Eq.  (32).  Collecting 

the  individual  Fourier  components,  one  obtains  four  separate 

systems  of  equations  for  the  constant  a i1?,  i  =  1,...4.  We  combine 
(i) 

the  constants  a^  ^  (for  fixed  values  of  i  and  k)  into  a  vector 
^  (whose  it,*"*1  component  is  a^^).  Then  these  problems  can  be 
written  in  the  form 


M(l)  a^l}  -  X^x)  0 


(35) 


L 


the  matrices  M ^ 1 ^  are  given  by 

,(i)  _  r>(i) 


M 


=  D'-1-'  +  (u2/4)M(l'1) 


(i)  , 


where  the  D  1  s  are  diagonal  matrices  with 


Dll)  =  ( 2n-2 ) 2 
nn 


d(2)  =  2 
U11  ^ 


Dnn}  =  ^2-^2 


-  '2"»2 


d(4)  =  o 
11  u 


D(4)  =  (2n-l) 2 
nn 


n  =  1,2,.. 


n  = 


n  = 


2,3,.. 

1,2,3, 


n  =  2,3, 


The  matrices  M 


(i,l) 


are  symmetric.  One  has 


=  2 1//2 

12 

=  1  , 

n,n+l  n+l,n 


n  =  2,3,.. 


(36) 
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* 


1 


=  1 

n,n+l  n+l,n 


n  1 f  2  f  •  »  9  f  l 

(36)  continued 


all  other  elements  of  the  matrices  are  zero. 

Let  A  ^ 1  ^  be  a  matrix  whose  k*1*"1  column  is  given  by  and 

let  be  a  diagonal  matrix  with  the  k^  element  X^*^  •  Then 

one  can  write 


M(i)A(i)  _  A(i)AU)  =  0 


I f  an  eigenvalue  routine  for  matrices  is  available,  (best  a  routine 

for  banded  matrices) ,  then  numerical  solutions  to  the  problem 

(Eq.  (35))  are  readily  obtained  by  truncating  the  infinite  matrix. 

In  the  present  context,  only  a  limited  number,  n,  of  eigenvalues 

and  eigenvectors,  a^1  is  used.  The  infinite  matrix  is  truncated 

to  an  m  by  m  matrix  where  m  is  chosen  large  enough  so  that  the 

coefficients  a  ^ (k  =  1, ?.,... n)  is  smaller  than  an  assigned  small 
ni  ^  k 

number  •_  ,  while  the  largest  coefficient  is  normalized  to  1,  say. 

In  our  computations,  carried  out  with  n  =  5  and  c  =  10-5,  it  was 
sufficient  to  take  m  =  9  if  p  •:  4  . 

The  Mathieu  functions  are  nearly  monotonic  (not  more  than 

one  minimum  of  the  absolute  value  and  no  maxima)  in  the  region 
2 

where  ..  -  (;:  /4)  cos  2n  •  0;  they  are  oscillatory  where  this 

2 

quantity  is  larger  than  zero.  If  u  is  not  small,  one  has  a  nearly 
monotonic  region  for  the  first  eigenvalues  in  the  vicinity  of  n  =  0 
and  n  =  r.  In  the  vicinity  of  n  =  n/2  and  n  =  3t;/2  the  functions 
are  always  oscillatory.  The  oscillatory  regions  grow  as  ?  is 
increased . 

For  the  separation  constants  obtained  from  the 

eigenvalue  problem  one  must  solve  Eq .  (31)  for  the  radial  Mathieu 

functions.  For  r  sufficiently  larae  these  functions  are  always 

2 

oscillatory,  because  contribution  (u  /2 )  cosh  (2r)  makes  the 
coefficient  of  f  positive.  The  particular  solutions  which 
represent  outgoing  waves  are  found  by  studying  the  asymptotic 
behavior  for  large  values  of  r.  It  is  practical  to  replace 
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by  u,  (Eq.  (23) ) ,  Then  one  obtains  from  Eq.  (31) . 


d2f1^l)  .  di'^l)  2 

_ +  l  — L  _  +  (if_ 

,  2  u  du  1  4 

du 


ai!X)  2  n 

k  p  lf(i) 

2  4 '  k 

u  4u 


(37) 


The  form  of  the  dominant  part  of  the  coefficient  of  f^'  suggests 
a  hypothesis 


f^X)  (u)  =  f^l)(u)exp(-i(u/2)u) 


(38) 


Outgoing  waves  are  obtained  by  choosing  here  the  negative  sign  in 
argument  of  the  exponential  function,  because  the  argument  of  the 
exponential  function  in  the  hypothesis  (3)  has  the  positive  sign. 

Then  one  obtains 


d2f 


(i) 


du 


df 


(i) 


ip) 


du 


1  2u 


(39) 


In  subsequent  formulae  the  indices  k  and  i  will  be  omitted.  A 
formal  development 


oo 


0 


c  u 
n 


-n 


(40) 


leads  to  the  recurrence  relation  (obtained  by  collecting  terms 

....  -(1/2)  -  (n-1 ) , 

with  power  u  ), 


ipn  cn  =  cn-1(A  -  (1/2  -  n)2)  -  (y2/4)cn_3  (41) 

According  to  the  hypothesis  (Eq.  (40))  all  coefficients  with 

negative  subscripts  n  are  zero.  The  last  equation  is  therefore 

-1/2 

satisfied  for  n  =  0.  (This  shows  that  the  factor  u  in  Eq.  (40) 

is  correct.)  The  coefficient  Cp  is  taken  to  be  1.  From  there  on, 
one  proceeds  to  larger  values  of  n.  The  series  is  semiconvergent . 

If  it  is  terminated  after  a  finite  number  of  terms,  then  the  error 
can  be  expected  to  be  of  the  order  of  the  first  neglected  term. 


This  determines  how  many  terms  and  which  minimum  value  of  u  one 
must  take  in  order  to  obtain  an  assigned  accuracy.  This  procedure 
works  well  provided  that  the  value  of  y  is  not  too  small. 

For  a  value  of  u  chosen  according  to  this  discussion  one 

computes  f  and  df/du.  With  these  initial  values  one  then  integrates 

Eq.  (37)  in  the  direction  toward  smaller  values  of  u.  For  larger 

values  of  u,  the  values  of  f  and  df/du  can  be  computed  directly 

from  the  asymptotic  representation  derived  here.  Notice  that  f, 

k 

and  df^/du  are  complex. 

Eq.  (37)  is  linear,  it  therefore  can  be  reduced  to  a  first 
order  differential  equation.  One  sets 

p(u)  =  (df/du) /f (u) 

or  (42) 

df/du  =  p (u)  f  (u) 

Hence 

d^f/du2  =  (dp/du)f(u)  +  p(u) (df/du) 

=  [(dp/du)  +  p2(u)]f(u) 


This  procedure  may  be  numerically  advantageous  especially  if  only 
(df/du) /f(u)  is  needed,  for  the  function  p  is  likely  to  be  smoother 
than  the  function  f.  But  this  is  merely  a  technicality. 
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In  the  present  problem  the  particular  solutions  in  terms 
of  Mathieu  functions  are  natural  because  a  series  in  these  particular 
solutions  (with  coefficients  determined  by  the  boundary  conditions 
at  the  plate)  converges  throughout  the  flow  field.  One  thus  obtains 
a  representation  of  the  flow  field  and  a  formulation  of  the  far 
field  conditions  with  which  other  approximations  can  be  compared. 

To  a  large  extent  this  observation  is  valid  also  for  problems 
in  which  the  equation  obtained  by  transforming  the  basic  equation 
differs  somewhat  from  the  Helmholtz  equation.  The  use  of  a 
representation  of  the  field  in  terms  of  Mathieu  functions  (n) 

need  not  be  restricted  to  theoretical  discussions.  For  practical 
applications  one  must,  of  course,  first  provide  representations 
of  the  functions  (n)  and  of  f^^  (u)  /df^j'du)  for  the  value 

of  u  at  which  the  far  field  conditions  are  applied.  This  can  be 
done  by  means  of  the  relations  derived  above. 

In  the  region  r  >  1  the  particular  solutions  ,  (Eq.  (30)) 

can  be  expressed  in  terms  of  the  particular  solution 
Eq.  (29).  In  principle  one  will  evaluate  a  particular  solution 

along  some  curve  i  =  const  (r  >  1),  and  then  represent  it 
in  terms  of  the  functions  x.p  by  means  of  a  Fourier  analysis  with 
respect  to  0 .  The  result  has  the  form 


;  i )  _  „  ( i ) . 

:  "  <LS’,kxJl 


Certain  coefficients  will  be  zero  because  of  symmetry  properties. 

The  coefficients  3  ,  are  constant.  The  values  of  the  constants 

(not  needed  in  the  future  discussions)  are  readily  obtained  if 
one  carries  out  the  above  matching  procedure  along  a  line  of  u  =  const 
in  the  limit  u  ►  -.  Consider,  for  instance  1  .  One  has, 

according  to  the  first  of  Eqs.  (34),  using  the  representation  of 
f,  given  by  Eqs.  (38)  and  (40) 

K 

=  g^1}  (n)  f,U)  (u) 


cos  (  (28.-2)  n)  ]  u  ^^exp  ( -i  (  p/2)  u) 
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we  write  Eq .  (4b)  in  a  form  where  the  notation  for  the  coefficient, 

6  is  analogous  to  that  for  the  coefficient  in  "a"  in  the  last 
equation 


B 


(1) 

l.k 


Hq  (Mr)  +  ^ 

l.J 

9=2 


,(D 

’e,k 


H 


(2) 

29-2 


(pr)cos  (  (2  9-2)  (•) 


One  has  for  u  large 


r  "e  u/2,  0  %  n 


ana 


l  7 ) 

«k  1  '  r  )  exp  [  i  (  r/2  )  (k  +  ]  (  2/ ■:  yr )  1/2exp  ( -  i  v.r  ) 

One  then  obtains 


6,'1.1  =  a.'1.' 


l,k  “l,k 


(i/2)  (hit)  exp  (-i  (tt/4) 


\  ( ^  )  __  (1)  1/0 
1 1  k  aJ. ,  k  f'^-/2)(i.  -’)  exp(-i 


(  tt/2  )  (  2  £  -  2  ■  Y 


>)] 


(46) 
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SECTION  V 

EVALUATION  OF  THE  FAR  FIELD  CONDITIONS 


Our  discussions  will  be  based  on  a  discretization  of  the 
differential  equation  in  the  u-n  system.  Accordingly,  one  must 
express  the  far  field  conditions  (Eq.  (28))  in  terms  of  these 
variables.  One  has 


=  gradcf  •  e  =  f  cos(n,x)  +  ;■  sin(n,x) 
n  n  x  v 

and 

ds  =  -dx  sin(n,x)  +  dy  cos(n,x) 

!  ds  =  e  dy  -  &■  dx  (47) 

n  x  2  v  v 


where  e  is  the  unit  vector  normal  to  the  contour  C;  (n,x)  denotes 
the  angle  between  the  normal  and  the  x-axis,  and  dx  and  dy  are 
the  components  of  a  line  element  along  the  curve  C,  which  is  now- 
given  by  an  ellipse  u  =  const.  Then,  from  Eq .  (19) 


dx  =  -Im(g ' ) dn 
dy  =  Re(g')dn 


y 


n 


(48) 


(49) 


Substituting  Eq.  (48)  and  (49)  into  Eq.  (47)  and  using  Eqs.  (20) 
and  (23) ,  one  finds  that  along  a  line  u  =  const 

$  ds  =  <J)  ^dn  =  U(|)  dn  (50) 

n  t,  u 


1 


Analogously, 


iij  ds  =  um  dn 
n  u 
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If  both  <p  and  w  are  considered  as  functions  of  u  and  r\ ,  then 
one  obtains  as  far  field  conditions 


w  d> )  d  n  =  o 

u 


(bl) 


Substituting  here  the  expressions  i p,  ,  Eq .  (30),  one  obtains 


fk<u>  /  ak(n)4>u(u,n)dn  -  fk(u)  /  gk  (n)  <t>  (u,  n)  dn  =  o 


where  the  integration  is  to  be  carried  out  along  the  line  u  =  const 
which  gives  the  contour  C.  The  values  of  fk(u),  fk(u)  and  the 
functions  gk(n)  are  considered  as  known. 

In  an  alternative  derivation  of  the  last  equation  one  would 
carry  out  a  Fourier  decomposition  of  <j>u(u,ri)  and  <j>(u,n)  (u  =  const) 
in  terms  of  the  system  of  orthogonal  functions  9k(n)  and  postulate 
that  each  of  the  components  represent  an  outgoing  wave.  This 
derivation  may  appear  simpler,  but  it  is  applicable  only  if  the 
contour  C  is  a  line  u  =  const  while  in  principle  the  formulation 
(28)  is  valid  for  any  contour. 

Further  formulae  are  needed,  if  one  chooses  for  the  test 
functions  w  in  Eq.  (28)  the  expressions  xm  defined  in  Eq .  (29) 

Xm  =  H^2) (pr)  cos(mO) 

We  assume  that  a  routine  for  generating  Bessel  functions  of  real 

argument  is  available.  Let  Y  (z)  denote  Bessel  functions  of  the 

m 

second  kind.  One  has 

Hm2’  *  Vz)  -  1  Ym(zl  (53> 
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We  note  the  following  relations 


H(2)  (z)  '  =  dH(2,/dz  =  1/2  H^2)  (z) 
m  in  x 


1/2  H<2>(z, 


Hm-l(z)  +  Hm+l(z)  =  (2m/z)HiJi2)  (z) 


In  evaluating 


f  y  d>  ds 

J  Am  yn 


one  uses,  again,  Eq .  (50)  to  express  cf>  ds.  The  radius  r  and  the 

angle  6,  which  are  needed  in  the  evaluation  of  Xm  are  expressed 
in  terms  of  n  (and  u)  by  Eqs.  (15)  and  (16),  the  trigonometric 
functions  of  0  ultimately  by  Eq.  (18).  It  is  best  first  to  compute 

C  =  £  +  in  =  log  u  +  in 

and  then  to  proceed  with  complex  arithmetic  (along  the  curve  C 
for  which  the  evaluation  is  carried  out,  one  has  u  =  const.) 

In  addition,  one  must  evaluate  (dxm/du)ds.  One  has 


dH(Vds 


37<Xm)r  d(-  -  —  (XJ  k  dr 


3  0  Am  r 


With  Eq .  (21)  one  then  obtains  for  a  line  u  =  const 


£-( ■  Ids  = 
dn  m 


rF  <Am)r  Re{(T}  +  3T{xm)  Im  ^  dn 


Hence,  after  substitution  of  the  specific  expressions  (29) 


■  Ids  =  {  r : ;  H  *  ^os  (m.O  )  Re  (— )  -  m  H  ^  2  ^  s  in  (m9  )  Im  (2--)  •  dn 
dn  m  m  g  m  g 


An  alternative,  more  symmetric  form,  is  obtained  by  using 
Eqs.  (54) ,  (55) ,  and  (22) 
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sin ( (m-1 )  ‘  ) 1m (q ' )  ] 


ak(Vds  = 


H  *  2j  ( \ir )  |cos  (  (m-1 )  •' )  Re  (q  ’ )  - 
m-i 


( 2 ) 

+  |  (  nr )  |  cos  (  (m+1 )  •  )  Re  (q  * }  +  sin((m+l)  )Im(g')l  d- 


{‘il) 


Now  the  integral  Eq.  (51),  with  u;  given  by  «  can  be  evaluated.  One 
needs  for  this  purpose  Eqs.  (50)  and  (56)  or  (57). 

A  formula  for  the  anproximate  far  field  conditions  of  Baydiss, 
Gunzberger  and  Turkel  is  obtained  by  the  following  considerations. 
The  fact  that  particular  solutions  representing  outgoing  waves  in 
terms  of  Mathieu  functions  have  the  asymptotic  form  for  large 
values  of  u 


y'k(u,n)  =  gk  ( n)  exp(-i  (|i/2)  u)  (u  1//2  +  c1  k  u-3^2  +  ..,) 

(see  Eqs.  (38)  and  (40))  imoiies  that  (  is  asymptotically  given 
by 

<f>  =  exp(-i  (ij/2)  u)  [u  1//2  h-jtn)  +  u  3//2  h2  ( '  )  +  ...1  (58) 


where  h^ ( n )  is  arbitrary  and  function  Iv  ,  i  1  are  expressed  in 
terms  of  functions  h^  and  their  derivatives  with  subscripts  j  •  i. 
The  operation 


applied  to  the  first  two  terms  in  Eq .  (58)  gives  zero  for  every 
value  of  n.  Writing  this  operation  in  detail,  one  obtains 


p  +4>  (—  +  ip) 
Kuu  U  U 


+ 


A  ill 

2  u 


0 


(59) 


This  equation  is  satisfied  for  any  function  <t>  which  has  the  form 

2  2 

of  the  first  two  terms  in  Eq.  (58).  The  term  1  <}>/3u  ,  which  does 
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not  fit  the  concept  of  a  boundary  condition  for  a  second  order 
partial  differential  equation  is  eliminated  by  means  of  the 
governing  Eq .  (27).  Then  one  obtains 


(60) 

This  condition  is  applied  along  a  line  u  =  const  for  0  £  n  5  2 tt  . 

The  presence  of  <j>  can  be  taken  as  an  indication  that  this 

nn 

formulation,  although  it  has  local  character,  is  an  approximation 
to  conditions  of  a  global  nature.  If  one  forms  analogous  conditions 
on  the  basis  of  additional  terms  in  Eq.  (58),  then  one  encounters 
higher  derivatives  with  respect  to  n. 
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SECTION  VI 

RELATION  BETWEEN  DIFFERENT  FAR  FIELD  CONDITIONS 


It  is  convenient  to  symbolize  the  function  f  (n),  f  (i),  .  (r  ) 

(n)»  \  ( n, )  ,  \  (n)  and  cos  mn  (or  sin  ran),  which  are  actually 

elements  of  a  function  space,  by  vectors.  The  integrations  needed 
in  the  far  field  conditions  are  then  considered  as  scalar  products 
which  will  be  written  as  products  of  a  row  and  a  column  matrix. 

One  of  the  far  field  conditions  then  appears  in  the  form 


[♦  ..  (U,r.)  ~  ['  m  (u,n)  ->]  {"  ‘  | 

.  I-  (u,  n) 1  -  u  U)  (u,  n) 

i 

u  f  j  L  ;  J 


o 


Here  one  must  substitute  for  u(n)  all  functions  of  either  the  set 
.,  (n)  or  the  set  >,  (n).  Combining  symbolically  these  row  vectors 
,  into  matrices,  one  then  obtains 


and 
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A 
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u  4  J 


r.  ...  -."i  p  , 

-Xk,u(u,n)  !;(u,-) 

■<-...  -+.J  L  j 
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In  all  of  these  equations  u  is  constant.  Because  of  the  special 
form  of  the  functions  4'^  (Eg.  (30)),  one  can  rewrite  Eg.  (61): 


D.  B  <J>  (u ,  n)  ’ 


D,  B 


1  ( u  ,  n  )  I 


=  0 


(b: 


where  and  D2  are  diagonal  matrices  in  which  the  elements  D 
and  D2  ^  are  respectively  given  by  f^(u)  and  f £ ( u 1  and  where 


B  = 


•  •  •  >“ 

♦-  gk(n)  ' 


Because  of  the  orthogonality  of  the  functions  g^fn)  one 
that  the  eigenvectors  p^fn)  of  the  eigenvalue  problem 


“  t 

-A 

Pk(n) 

~  X.D-B 
k  2 

Pk(n) 

m  't 

.  1 

are  given  by 


pk(n)  =  gk{,i) 


and  that  the  eigenvalues  are 


k 


Dj_,kk 

D2,kk 


fk(u) 

f^lul 


The  relation  between  the  formulations  Kq .  (63)  and  (62)  is 

established  on  the  basis  of  Kg.  (4  )  which  is  now  written  in 
form 


where  ^  is  a  vectoi  in  the  ordinary  sense, 
component  is  given  by  .  Hence 

'  r  K 


for  which  the 


th 


,  kk 


(64) 


finds 


(65) 


(66 ) 


(67) 


he 


where  r  denotes  a  constant  matrix  in  the  usual  sense  for  which 
the  element  with  subscripts  k  and  •  is  given  by  „.  One-  rccoun  i  zes , 
by  substituting  the  last  relation  into  Eq .  (61) ,  that  Eq .  (62)  arises 

by  premultiplying  Eq.  (61)  by  the  matrix  ((']-  The  matrix  Irl  is  well 
behaved.  Tne  relative  eigenvectors  and  eigenvalues  of  the  two 
matrices  occurring  in  Eqs.  (61)  and  (62)  are  therefore  the  same. 

After  a  discretization  this  result  will  hold  only  approximately, 
the  deviations  from  the  ideal  case  will  serve  to  characterize  the 
effect  of  the  discretization. 

In  discussing  the  approximate  far  field  conditions  of  Bayliss, 
Gunzburger,  and  Turkel,  Eq .  (60),  it  suffices  if  one  considers  one  com¬ 

ponent  of  the  development  of  i(u,n)  and  if  (u.t<)  along  a  line  u  =  const 
with  respect  to  the  functions  g^(n).  We  assume  accordingly  that 


<$»  <u,  n)  =  qk  (u)  gk  (  n) 


then 


(69; 


d>u ( u , n )  =  q£(u)  gk(n) 


The  function  qk(u)  replaces  the  function  f^tu)  in  Eq .  (30).  Then 
one  obtains  from  the  approximate  far  field  conditions  Eq .  (60) 


da 


;k 

du 
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+  in) 
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1_ 

'  2  ,2 

U  dr- 


( COS  2 r 
4u  2u 


2  .  2 

T  ~  7~4}ak] 

4  u 


and  by  substitutincj  Eq .  (32) 


^k 

du 


+ 


3  i_a 
2  u 


0 


(70) 


The  same  relation  is  obtaincci  if  one  assumes  that  a,  has  the  form 

•k 

a,  ( u )  =  (  ,  u  *  2  +  u  1  2 )  exo  ( -i  ( .  ,/2 )  u) 
k  i  2 


where  and  a are  arbitrary  constants.  Any  expression  of  this 
form  satisfies  the  second  order  differential  equation 


(u) 


0 


In  the  derivation  of  Eq.  (60) ,  the  condition  has  been  imposed 

that  <t>  (u , n )  (here  given  by  Eq.  (69))  satisfy  the  partial 

differential  equation  (27)  at  the  value  of  u  which  gives  the 

curve  C.  This  leads  to  Eq.  (37)  with  fk(u)  replaced  by 

q,(u)  .  Eq.  (70)  is  then  obtained  by  eliminating  from  the  last 
k  2  ? 

expression  d  qk(u)/du“  by  means  of  Eq.  (37). 

The  eigenfunctions  pertaining  to  the  approximate  conditions 
Eq .  (60)  are  again  the  f unctions  g^  (  q)  >  but  the  eigenvalues 
fk (u) /f£ (u)  are  replaced  by  the  expression  obtained  from  Eq.  (70). 
This  equation  arises  by  assuming  that  the  function  gk(u)  which 
replaces  fk(u)  has  the  correct  asymptotic  form,  and  satisfies 
the  differential  equation  for  f  for  the  chosen  finite  value  of  u. 


SECTION  VI 1 
DISCRETIZATION 


So  far  the  far  field  conditions  have  been  formulated  as  if 
the  functions  f  involved  can  be  represented  exactly,  and,  at 
Least  in  principle,  it  has  been  assumed  that  the  set  of  test 
functions  to  be  used  in  the  far  field  conditions  is  infinite. 
Actually,  the  functions  <j>  must  be  approximated  by  means  of  a 
finite  number  of  expressions  and  the  set  of  test  functions  for 
the  far  field  conditions  can  only  be  finite. 

An  ideal  case  arises  if  <f>  satisfies  the  Helmholtz  equation 
even  in  the  computed  part  of  the  flow  field.  Of  course,  one 
would  then  proceed  analytically,  but  this  case  serves  well  in 
comparing  different  formulations.  If  one  represents  the  solution 
:  by  a  superposition  of  particular  solutions  ; ^  (defined  in 
Lq .  (30)),  the  Helmholtz  equation  and  the  far  field  conditions 

are  satisfied.  To  satisfy  the  boundary  conditions  at  the  plate 
one  must  use  an  infinite  series  in  these  particular  solutions. 

This  series  converqes  throughout  the  flow  field.  In  practice,  one 
will  truncate  the  series.  The  only  error  which  arises  here  is 
caused  by  the  approximation  to  the  boundary  conditions  at  the 
plate . 

In  cases  of  practical  interest,  the  equation  obtained  by 
transforming  che  original  problem  differs  from  the  Helmholtz 
equation,  particularly  in  the  vicinity  of  the  plate.  (It  is  of 
course  assumed  that  at  great  distances  the  transformed  equation 
is  very  close  to  the  Helmholtz  equation).  The  transformation  of 
the  flow  field  to  a  half  infinite  strip  in  the  ’  ,  r,  plane  may 
nevertheless  be  practical,  as  the  transformation  provides  a 
suitable  mesh  in  the  physical  plane.  The  fact  that  in  this  report 
the  analysis  is  carried  out  in  the  r  ,  t  (or  rather  in  the  u,n) 
system  ought  to  be  considered  as  incidental;  it  facilitates  the 
integrations  necessary  in  the  formulation  of  the  far  field 
cond i t ions . 
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1 


In  the  u,n~plane,  <J>  is  periodic  along  lines  u  =  const.  For 
a  function  which,  in  addition,  is  analytic  along  such  lines,  a 
representation  in  terms  of  analytic  periodic  functions  has  excellent 
convergence  properties.  (The  coefficients  a^(sav)  of  a  Fourier 
series  decrease  for  sufficiently  large  n  faster  than  any  negative 
power  of  n.)  One  such  representation  is  by  means  of  a  truncated 
series  in  Mathieu  functions 


where,  so  far,  the  coefficients  q^  are  unknown  functions  of  u. 

If  the  problem  is  governed  by  the  Helmholtz  equation,  then  this 
hypothesis  leads  back  to  a  representation  of  <J>  by  a  superposition 
of  functions  iji,  .  The  only  errors  arise  because  general  boundary 
conditions  at  the  plate  cannot  be  satisfied  exactly  by  a  truncated 
series . 

For  a  more  general  differential  equation  one  is  led  to  a 
system  of  second  order  ordinary  differential  ecuations  for  the 
functions  q^(u).  At  large  distances  the  functions  qk(u)  will  have 
the  same  behavior  as  the  functions  f^(u),  which  is  determined 
by  Eq.  (37) .  For  u  large, one  has  approximately 

fk  =  u"1/2  exp(-i (u/2) u) 

This  means  that  the  relative  maanitude  of  these  functions  for 
different  values  of  k  remains  the  same  (unless  k  is  very  large, 
then  the  asymptotic  development  is  not  applicable).  Therefore, 
the  number  of  functions  gk(n)  in  the  representation  of  1 
cannot  be  reduced  as  one  proceeds  from  intermediate  to  large 
values  of  u.  Such  a  reduction  would  be  possible  for  the  Laplace 
equation.  More  functions  gk(n)  may  be  needed  close  to  the  plate 
in  order  to  express  local  effects  of  the  boundary  conditions. 

To  reduce  from  a  certain  value  of  u  on  the  number 
of  functions  g^ ,  one  will  set  functions  g^tu)  for  the 
expressions  that  are  to  be  disregarded,  equal  to  zero.  This 
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implies  that  at  this  value  of  u  one  introduces  the  condition  =  0 
for  those  values  of  q  .  The  functions  q,  (u)  resemble  the  functions 

K  2  2 

f^  (u)  and  those  functions  are  nearly  monotonic,  for  n  ;i  u  /4  <  >  , 

for  larger  values  of  u,  they  are  oscillatory.  If  one  imposes  the 
condition  q^(u)  =  0  for  a  value  of  u  where  the  function  f^  (u)  is 
oscillatory,  then  there  is  a  possibility  that  large  errors  are 
introduced  because  of  a  resonance  phenomenon.  To  avoid  this,  one 
will  reduce  the  number  of  functions  g^fn)  used  for  the  representation 
of  <J>  only  at  values  of  u  for  which,  for  the  value  of  k  in  question, 
f^  is  nonoscillatory . 

If  <j)  is  represented  at  each  value  of  u  in  terms  of  Mathieu 
functions  g^(n),  then  one  must  provide  a  representation  for  these 
functions,  obtained  for  instance  by  the  procedure  shown  in  Section 
IV.  This  may  well  be  a  worthwhile  approach.  In  this  case,  no 
approximations  ar  :  encountered  in  the  far  field  conditions. 

The  periodicity  of  1  with  respect  to  n  is  alternatively 
expressed  by  the  use  of  a  truncated  Fourier  series 


<Mu,n)  ~  >  q,  (u)  hk(>:) 


where 


hQ(n) 


2_1/2 


(72) 


h0k(n)  -  cos(kr-),  hi  2  ^  _  i  {  r, )  =  sin(kn);  k  =  l,...2 

In  the  hypothesis  (72)  a  certain  minimum  number  of  functions  g^O;) 
are  needed  in  order  to  give  a  sufficiently  accurate  representation 
of  1  for  large  values  of  u.  We  have  also  shown  that  the  functions 
g  (n )  can  be  represented  by  a  Fourier  series  in  .  In  this  sense 
the  hypothesis  Eqs.  (71)  and  (72)  dire  related  to  each  other.  All 
Fourier  components  (implicitly)  present  in  Eq .  (71)  must  also  be 

present  in  Eq .  (72)  in  order  to  give  a  desired  accuracy  for  large 

values  of  u.  This  means  that  the  maximum  value  of  k  occurring  in 
Eq.  (72)  must  exceed  that  occurring  in  Eq .  (71) .  The  number  of 
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additional  terms  depends  upon  the  value  of  p ,  it  increases  with 
[i .  For  the  detailed  discussion  of  this  report  we  have  used  the 
function  h  defined  in  Eq .  (72).  Sometimes,  however,  one  might 

find  it  desirable  to  use  different  functions  for  h^Oi)  .  The  good 
convergence  properties  of  the  Fourier  series  hold  only  if  the 
functions  to  be  represented  are  analytic  along  a  line  u  =  const. 

If  they  have  discontinuities  (for  instance  in  higher  derivative) 
or  even  if  they  are  analytic  and  strongly  peaked  (that  is,  if  in 
the  complex  n-plane  poles  lie  close  to  the  real  n  axis)  then  the 
convergence  will  be  poor.  For  subsonic  flows  the  governing 
differential  equation  has  analytic  coefficients  and  f  will  be 
analytic  except  at  the  plate  where  the  boundary  conditions  may 
introduce  singularities.  For  flows  with  an  embedded  supersonic 
region  even  the  coefficients  of  the  differential  equation  need 
not  be  analytic  (especially  toward  the  end  of  the  supersonic 
region) .  Under  these  circumstances  one  will  prefer  representations 
which  are  more  suitable  to  reflect  local  properties  of  $ ,  for 
instance  finite  difference  or  finite  element  approaches.  This 
expresses  itself  in  the  choice  of  the  functions  h^fn)  •  For  the 
value  of  u,  at  which  the  far  field  conditions  are  applied,  one 
deals,  of  course,  with  analytic  functions.  Adopting  the  notation 
of  the  preceedina  section,  we  write  the  expression  Eq .  (72)  in 

the  form 

r f  n 

where  the  number  of  functions  hr(n)  and  of  components  of  the 

A 

vector  q(u)  is  finite.  The  expression  Eq.  (71)  arises  by  replacing 
in  this  equation  hk(n)  by  gk(n). 

Let  us  first  study  Eq.  (61) ;  in  which  the  far  field 
conditions  are  expressed  by  means  of  Mathieu  function.  They  are 
written  in  the  form  of  Eq .  (63). 
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The  number  of  equations  expressed  by  this  formula  equals  the-  number  of 
rows  in  the  matrices  D^B  and  D2B.  This  number  must  be  chosen  equal 
to  the  number  of  components  in  the  vector  q(u).  Let  this  number  be 
given  by  n^.  We  substitute  info  this  equation  tsc  expression  (64)  foj 
B .  Let 
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r  "j 

r  *  4  ♦ 

1  •  » 
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i 
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Let  Dj  and  D2  be  the  matrices  and  D-,  and  D9  truncated  in  the  suim 
manner.  Then  one  has 


D1  B1 


q' (u) i 


D2  q ( u ) 


=  0 


(73) 


In  the  preceding  section  we  had  characteri zed  the  far  field 
conditions  by  the  relative  eigenvectors  and  eigenvalues  of  the  two 
governing  matrices  (here  D-^  B-^  and  D~,  B^)  .  Accordingly,  we  consider 
the  problem 

~  vf  °2  "  0  (76) 

where  v0  is  the  eigenvalue  and  Vf  the  ?fc  eiaenvector.  Obviously, 

V,  =  B. -1  e„  (77) 

y.  1  k 

th 

where  e.  is  the  unit  vector  in  the  direction  of  the  •  component  o! 

V  (i.e.,  the  ktIri  component  of  e,  is  given  by  f ,  ,).  Then 

*  ‘  ,  K 

v?  ®l,  9  9^2 ,  9.9  (  /  8' 
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This  result  holds 
I  f  )  =  dk  (’!)  , 

t  hi'  m.i  t  r  i  x  B  j  is 
i  re  norma  1 i zed )  . 


whether  is  an  infinite  or  a  finite  matrix, 
then  because  of  the  orthogonality  of  the  q^(n) 
the  unit  matrix  (provided  that  the  functions  cj^,  ( r, ) 
Then 


v 

9. 


The  function  f(u,r))  pertaining  to  this  eigenvector  is  given  by 


t  t  1 1 

.  .  .  qk(n)  .  .  .. 

}  4  4-  ! 


+ 

g^n) 
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This  is  the  result  obtained  in  Eq .  (66).  It  holds  whether  the  matrix 

B  is  truncated  or  not.  In  this  case  no  error  will  arise  in  the 
far  field  conditions. 

Assume  next  that  the  functions  h,,  (n)  are  the  terms  of  a 
••'ourier  series.  Let  be  the  matrix  B  ^  without  truncation.  Then 
■■me  finds  in  the  row  of  the  matrix  the  Fourier  coefficients 

of  u..  The  eigenfunctions  are  given,  as  before,  by  the 

1  - 1  . 

column  of  .  For  the  nontruncated  matrix 


B 


1 


-1 


T 


(where  T  stands  for  transpose).  This  is  seen  as  follows.  The 
mat  r  i  x 

t  t  ? 

.  .  .  .  h .  ( n )  •  -  : 

l 

i  i  } 


maps  an  element  of  the  vector  space  (in  which  V  lies)  into  a 
function  space  spanned  by  the  individual  functions  h.  (•■)  . 
Specifically,  it  generates  the  functions  which  has  the  components 
>■•:  V  as  coefficients  of  the  Fourier  decomposition.  On  the  other  hand 
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yields  the  coefficients  of  the  Fourier  decomposition  of  p(n) 
because  the  h^(n)  are  assumed  to  belonq  to  an  orthogonal  set 
It  follows  from  this  description  that 


p . H 
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and  on  the  other  hand  that 
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The  same  holds,  of  course,  for  any  other  orthonormal  set  q ,  (n) 

K 

It  then  follows,  because  of  Ea .  (74),  that 


B1  B1  = 


h,  (".) 


'  I 


qk(!) 


•h.  (’  ) 


is  the  .identity  operator.  Hence,  the  results  announced  above 


T  -1 

Bi  =  B1 

.  th 
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We  mentioned  that  the  j  1  row  of  gives  the  components  of  a 

decomposition  of  g.  in  terms  of  the  orthogonal  set  h(r).  According 

1  *■  ^  p 

to  Ea .  (77) ,  the  eigenvector  V-  is  given  by  the  j  column  of 

-IT  ^ 

B  ^  -  B  ^  .  This  means  that  for  nontruncated  matrices 
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*  i 

b  i 
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Tn  the  representation  of  the  flow  field  discussed  first, 

: (u,- )  is  expressed  for  each  value  of  u  as  a  linear  combine1 ion 
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of  Mathieu  functions  g£  (n)  .  The  eigenvector  V{  is  then  given 
by  the  unit  vector  e£  and  the  function  represented  by  it  is 
immediately  given  by  g?(n).  If  4>(u,n)  is  represented  for  each 

—y 

value  of  u  by  a  Fourier  series,  then  the  eigenvectors  Vj  has  as 
components  the  Fourier  coefficients  of  the  function  gj(n).  It 
follows  that  for  far  field  conditions  in  the  form  (63)  (i.e., 

far  field  conditions  which  use  Mathieu  functions  as  test 
functions)  the  eigenvectors  obtained  with  these  two  choices 
of  representations  for  0(u,n)  represent  the  same  function 
q,(v).  This  holds  only  for  nontruncated  matrices,  for  Eq .  (79) 

A. 

is  subject  to  this  restriction. 

T 

Eq.  (79)  is  based  upon  the  observation  that  B1  B1  is 

the  identity  operator.  For  completeness,  the  meaning  of  the 
—  T  — 

operator  B^  B^  (where  B^  denotes  the  truncated  matrix)  is 
explored  in  the  Appendix. 


In  practical  computations  one  need  not  evaluate  the 
eigenvectors  .  They  are  introduced  here  to  describe  the 
working  of  the  far  field  conditions.  One  has  according  to 
Eq.  (77) 


(80) 


Let  g?  be  the  function  that  is  represented  by  . 


g£(n) 


!  +  t  t]  j  ^ 

'  hk(n) 

1  4-  1  ; 


-> 


The  function  <j)(u,n)  can  be  expressed  as  a  superposition  of  these 
eigenfunctions 


Mu,n)  -  Iq£(u)  g£  (n)  ;  <f>u(u,n)  =  £q£(u)g£(n) 

Then  one  has  for  the  value  of  u  which  gives  the  curve  C 
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q ,  ( U ) /q ’  ( u ) 


=  v> 


f 


where  accordinq  to  Eqs  .  (78}  and  (67) 


\  ) 

'  £ 


?  f //D2 


fk(u)/f'k(u) 


(81  j 


For  truncated  matrices  the  functions  q^(n)  are  in  uencral 
not  identical  with  the  functions  q  (?■)  for  nont  runcutod  rntn-'c:;. 

K 

For  the  practical  application  of  the  far  field  condition.'. 

(75)  one  must  form  the  matrix  B1 .  (In  Eq .  (80)  its  inv-rso 

is  encountered.)  We  observed  above  that  the  row  of  B 

consists  of  the  Fourier  coefficients  of  the  function  >.r  .  (•  )  . 

1 

In  this  set  of  coefficients  is  truncated. 

The  coefficients  of  a  Fourier  development  of  the  Minot j 'ns 

q_j(r:)  have  been  used  previously  to  provide  ■  represent  ion  of 

these  Mathieu  functions.  They  arise  by  sol vi no  the  o iuenva 1 ue 

/  '  ) 

problem  Eq .  (35)  for  the  infinite  matrices  M'  .  In  nrac*  i  c*  ■ 

these  matrices  are  truncated  but  only  to  the  extent  Mi.it  a]) 
Fourier  coefficients  of  those  functions  n_.(>)  which  are  us*.  •! 
in  the  computations  are  found  with  a  prescribed  accuracy. 
Actually,  a  more  drastic  truncation  is  admissible  in  the  :  rts-nt 
context.  We  observed  above  that  the  number  of  terms  occur )- 1  nu 
in  the  sum  of  En .  (72)  is  larqer  if  one  represents  q  i  ,v  a 

Fourier  development  than  if  one  uses  Mathieu  functions  )  . 

The  last  few  eigenvectors  are  inaccurate  in  any  case  because 
of  the  truncation  of  B  .  This  can  be  tolerated  because  the 
number  of  functions  h^fn)  (of  Fourier  components)  is  supposed 
to  be  largo  enough,  so  that  in  the  (conceptual)  decomposition 
of  !■  the  last  few  eigenvectors  occur  only  with  very  small 
coef  ficients . 

The  accuracy  of  the  formulation  of  the  far  field  conditions 
need  not  exceed  the  accuracy  with  which  the  partial  differential 
equation  is  satisfied.  If  one  approximates  the  function  ‘ (up  ) 
by  a  truncated  series  of  the  form  (71)  then  the  resulting  system 
of  ordinary  differential  equations  is  simply  qiven  by  Eq .  (3.1) 

with  f^  replaced  by  . 
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If  one  represents  ;  ( u ,  •  )  in  the  form  (72)  (in  which  the 
:  ;:net  i  -  ms  h,  (•  )  arc  the  terras  of  a  Fourier  series)  then  one 
turns  hy  s’U'St  i  t  ntion  into  t  hi'  partial  differential  Eg .  (27) 

-I  *  u  ^q*  +  (..“/4)  (1+u  ^ )  q  -  u  ^Mq  =0  (82) 

*  t  h 

here  n(u)  is  a  vector  valued  function  whose  k  component  is 

uiven  by  the  function  q,  (u)  occurring  in  Eq .  (72).  The  matrix 

M  is  uiven  by  the  matrices  M  defined  in  Eq .  (36).  Some 

modification  of  the  indices  is  needed  because  the  Eqs .  (36)  take 

the  symmetries  of  the  problem  into  account. 

Fq .  '82)  holds,  in  the  first  place,  for  the  infinite  system. 

But,  because  we  use  in  Eq .  (72)  a  Fourier  series  truncated  to  n^ 

terms,  one  must  truncate  also  the  matrix  M  to  a  size  by  n^ . 

(I’)i  is  was  the  form  in  which  the  problem  has  originally  been 

or on rammed . ) 

The  matrix  encountered  in  the  far  field  conditions  (when 

on>-  uses  Mathieu  functions  as  test  functions  and  a  truncated 

Fourier  development  for  the  representation  of  j>)  ,  contains  in 
t  h 

its  j  row  the  truncated  set  of  Fourier  coefficients  for  the 

Mathieu  function  g . (n) .  They  are  computed  by  solving  the 

3  ( i  ) 

eiqenvalue  problem  for  the  infinite  matrices  M  ,  Eqs.  (36). 

This  has  been  stated  before. 

It  is  consistent  with  the  approximation  used  in  solving 

the  partial  differential  equation  if  one  approximates  the  rows 

of  the  matrix  B  by  the  components  of  the  eigenvectors  of  the 

( i  ) 

matrices  M  truncated  to  the  dimension  n^  by  n^. 

The  eigenvalue  problem  Eq .  (35)  yields  also  the  eigenvalues 

,  which  are  used  in  Eq .  (31)  to  determine  f,/f’  (Eq.  81)  for 

.<  K  K 

the  value  of  u  which  represents  the  outer  boundary  C.  The 

detailed  procedure  is  described  in  conjunction  with  Eqs.  (37) 

through  (41).  If  one  truncates  the  matrices  to  n^  by  n^ 

matrices,  then  there  will  be  an  error  in  the  later  eigenvalues 

an*1  in  the  values  f  /f'.  Such  an  error  is  admissible. 

^  K  K 
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The  error  encountered  in  this  procedure  can  he  v i run] i /ed 
if  one  computes  the  anqle  between  the  approximate  and  the  exact 
eigenvectors  g  (ri)  and  the  relative  error  in  the  ei  o.  n-M  ]  ues . 

K 

The  functions  g  (')  are  given  say  their  Fourier  eoet  !  i  c  1  •  ,-jg  ... , 
here  combined  into  vectors.  Consider  two  vectors  Vj  and  V 
with  components  and  b..  The  scalar  product  is  defined  by 

[Vj,  V21  -  a.  b.  (83) 

i 

The  angle  is  then  defined  by 


Ir.  the  present  case  the  vector  V_  which  represents  q.  (•  ) 

Z  '  K 

has  only  components.  Accordingly,  one  has  a.  -  0,  i  ■  iij  • 

■  ► 

The  vector  which  represent  g^fn)  has  an  infinite  number  of 
components,  but  the  Fourier  series  converges  very  we]].  The 
relative  error  in  the  eigenvalues  is  given  by 


The  approximate  far  field  conditions  of  Bayliss,  Gunzburqer 
and  Turkel  Eq .  (60)  arise  from  the  assumption  that  1  has  the  form 

of  Eq .  (58)  .  This  leads  to  Eq .  (59)  .  In  this  equation  is 

eliminated  by  means  of  the  original  partial  differential  equation 
Of  course,  the  requirement  that  the  partial  differential  equation 
be  satisfied  is  needed  only  for  the  value  of  u  which  corresponds 
to  the  outer  boundary  C.  The  approximation  to  the  condition 
(60)  for  cases  where  the  potential  is  approximated  by  the 
truncated  Fourier  series  Eq .  (72)  is  best  found  by  retracing  this 

derivation.  Eq .  (59)  then  assumes  the  form 


37 


I"  t.V  ( .hi'  1  >  1  .  )  f.'j  (-..J/4+  (  3/4  )  u  2*(3/2)i,.u  1  )•-•(),  03  5) 

Ft  til;;  n>is  q"  is  <•]  imin.iti'd  by  means  of  F.'<  i .  { 82 )  ,  which 

is  n.  .'Mj.'i  ess  i  on  of  t  ho  partial  differential  equation.  One 

obt a i ns 


ql  U>u'l  +  i.  )  *  i!-{  2  2)  (l*u  4  /  2  )  +  (  3  /  4  )  u  2+(3/2)i,.u  1  ]  +u  2Mq  =  0 

where  .•-us!  sulist  ltuto  flu'  value  of  u  for  the  contour  C.  This 

is  t  he  ' ' 1  it.  of  the  Bayliss,  Gunzburqer,  Turkel  condition  as  it 
is  list'd  in  tin-  practical  work. 

In  tins  equation  the  truncated  matrix  M  appears  aqnin.  For 
v’  mu  lysis  of  the  effect  of  this  boundary  condition  (in  eoniunct ion 
with  a  represent  a  t ion  of  t  ho  flow  field  by  a  truncated  Fourier 
series  in  • j we  assume  that  q  is  represented  by  a  linear  combination 
•  >f  the  i  t.  rtinent  oi convectors  V,  (which  in  turn  qivo  representations 
for  the  eigenfunctions  q  (>.)  .  bet  the  coefficient  do  q.  (u)  . 

Becuisi  of  the  definition  of  V.,  one  has 

M  V?  =  >f  Vf 

where  V.  denotes  the  truncated  matrix  M  and  >  f  is  the  eiqon value 
which  telonas  to  a  (’ ).  Therefore  from  Eq .  (85) 

o'  (  2  u  ~  1  m  .  )  to  [-  (..2/2)  ( 1  +u-4/2  )  +  (  ?/4)u_2  +  (3/2)  i  uu"  1  +  1  f  u~2 )  -0  (86) 

The  application  of  the  Bayliss,  Gunzburger,  Turkel  condition 
thus  has  the  following  effect.  For  the  value  of  u  con espond inq 
to  the  curve  C  one  decomposes  ;  ( u ,  n )  and  :  (u ,  ■  )  in  terms  of 
the  approximate  e igcnf unct ions  q  (n)  and  then  postulates  that 
the  ratio  of  the  coefficients  q'{  (u)  /qf  (u)  assumes  the  value 
computed  from  Eq .  (86).  The  error  encountered  in  this  procedure 

is  due  to  the  deviation  of  the  approximate  ciqen funct ions  a„(>) 
r  r '  >  m  the  exact  e icon funct ions  qf(n)  and  the  error  in  the  eigenvalues 
r  (u)  'f'  (u).  This  error  is  characterized  by 
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|  ff  (u)/ff  (u)  -  qv  (u)/qf  (u)  1 
!  f  .  (u) /f  (u)  ] 

X  X. 


(87) 


This  expression  must  be  evaluated  for  the  values  of  fc  used  in 
the  representation  of  the  solution  and  for  the  value  of  u  which 
gives  the  boundary  C. 

Finally,  we  study  the  effect  of  a  truncation  if  one  uses 
as  test  functions  the  expressions  <  defined  in  Eq .  (29).  They  are 

a  product  of  Hankel  functions  in  r  and  trigonometric  functions 
in  9  . 


It  was  shown  in  Section  VI  that  the  formulation  of  the  far 
field  conditions  in  terms  of  test  functions  originates  from  a 
formulation  in  terms  of  test  functions  b  by  premultiplication 
with  a  matrix  (Eqs.  (61),  (62),  and  (68)).  We  found,  furthermore, 

that  if  one  equation  of  the  system  Eq .  (62)  is  satisfied  for 

some  value  of  u,  it  will  be  satisfied  for  any  other  value  of  u. 


Nevertheless,  one  may  have  reservations  about  applying 
far  field  conditions  in  terms  of  functions  v  if  u  is  small.  Along 
each  line  u  =  const  the  function  x(u,n)  is  analytic,  and,  therefore, 
a  representation  in  terms  of  a  Fourier  series  in  n  is  possible. 

The  coefficient  of  Xj(u,n)  appear  in  the  j  ^  row  of  the  matrix 


-o 

1  t  t 

: . .  x .  ( u ,  n ) . . : ! 

h  .  (u) 

•  3 

3 

U  "J 

!  4-  , 

But  the  functions  x  are  singular  at  r  =  0;  one  has  Xm  ^  r  m  cos  m  . 
Along  a  line  u  =  const  one  has  small  values  of  r  for  n  =  r/2  and 
9  =  3t/2,  and  there  the  functions  x  have  peaks,  particularly 
pronounced  for  values  of  u  close  to  1.  Then  a  considerable 
number  of  Fourier  coefficients  will  be  needed  to  give  an 
acceptable  approximation  to  x •  Because  of  the  truncation  of  the 
Fourier  series  to  terms  the  coefficient  s  which  are  not  negligible 
will  be  disregarded. 
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Another  argument  seems  to  support  these  reservations.  The 
representation  of  a  test  function  in  terms  of  function  x  converges 
only  if  r 

For 


reached  at  the  y  axis.  There  r  =  y  =  sinh  f  =  5  (u  -  ^) 


1.  The  smallest  value  of  r  along  a  line  u  =  const  is 

1 
2 

u  •  1  +  /2,  there  exist  points  along  a  line  u  =  const  where  the 
development  of  the  functions  in  terms  of  Xj  does  not  converge. 
(The  reason  is  again  the  behavior  of  the  functions  Xy)  One  may 
then  have  doubt  whether  the  formulation  of  the  far  field  conditions 
in  terms  of  the  test  function  can  replace  the  formulation  in 
terms  of  test  functions  ^ •  Of  course,  the  underlyino  theory 
shows  that  the  functions  are  legitimate  test  functions. 


Substituting  the  hypothesis  (72)  (written  in  the  form  of 
Eg.  (73))  into  Eq.  (62),  one  obtains,  if  one  uses  nontruncated 
expressions 

!  I  L.  I 

0 


B2(u) 


■  + 

» 

,  q  ( u )! 


B3  (u) 


I 


q  (u) 


88) 


where 


B.. 


B 


3 


*  •  • 

*  •  xk  (u,n)->  ,  h  (-.) 

*  < 

t 

— 

A  A 

<  • 

\ 

j 

.  .Mr,) . . . : 

_ 

■ 

’ 

It  follows  from  Eq .  (68)  and  the  definition  (30)  for  1 p  that 

Eci .  (75)  arises  from  the  last  equation  by  premultiplication  with 

t  he  matrix  Its]  .  We  had  found  that  the  eigenvector  of  the 

matrix  DjB  relative  to  the  matrix  D2B  is  given  by  a  vector  whose 
components  are  the  coefficients  of  the  development  of  g.  with 

..  y  1 

respect  to  the  system  of  functions  fw  ,  and  that  the  j  eigen¬ 
value  is  given  by  f-(u)/f.(u).  The  eigenvalues  and  eigenvectors 

■J  J 

remain  unchanqed  by  the  premultiplication  by  [(B)]  .  These  are, 

therefore,  also  the  eigenvalues  and  eigenvectors  of  the  matrix  B2 
relative  to  B ^ .  Using  the  notation  of  Eq .  (35)  (but  without  the 

superscript  i,  which  is  introduced  to  distinguish  between  cases 
of  different  symmetries),  one,  therefore,  has 
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f  k  •  (  u) 

D,  B-  V  -  v  D„  B.  V,  -  0 

Ilk  (u)  2  ]  k 


and,  consequently,  also 


V(u) 

B2  Vk  fk(uT  B3  Vk  ° 


We  assert  that  for  the  first  few  eigenvectors  this  equ.it  ioi  i 
rather  well  satisfied,  even  if  the  matrices  B2  and  B,  .md 

•  -V 

vectors  V,  are  truncated.  This  assertion  is  based  on  the  act 
k 

that  the  truncation  of  the  first  few  eigenvectors  changes  tie 

eigenvectors  only  by  a  very  small  amount,  because  tho  comp'  n>  nt  h 

lost  by  the  truncation  are  close  to  zero.  The  trunoh  ion  of  B 

and  B ^  to  a  finite  number  n^  of  rows  means  simply,  that  one 

considers  only  the  first  n^  equations,  rather  than  all  of  them. 

The  truncation  to  a  finite  number  n^  of  columns  mi aht  .q.pour 

more  serious  because  the  elements  that  are  omitted  not 

small.  However,  in  the  original  infinite  matrix  they  ito 

multiplied  by  the  components  of  the  eigenvector  wh  i  •  -t .  ttn  sm  ill, 

and  which  are  now  eliminated  because  of  the  truncation. 

Incidentally,  the  eigenvectors  with  the  largest  ei  uenv.i  1  u< 

will  satisfy  the  truncated  equations  only  poorly,  because  1 1 . . •  •  i 

the  components  lost  by  truncation  are  lame.  We  emu),  as  i  :n  b  a!  o--, 

that  the  number  n^,  which  determines  the  size  of  the  mat!  i  .x ,  rust 

be  larger  than  the  number  of  eigenfunctions  which  are  important 

to  represent  the  solution  because  n ^  represents  the  functions. 

g^Oi),  which  are  important  for  the  solution  with  sufficient 

accuracy.  In  other  words,  for  n  '  n^,  the  components  V.  ^ 

»•  * 
with  n  '*  n^  of  the  vector  must  be  negligible  for  those 
*■ 

vectors  which  are  important  in  the  solution  !.  We,  therefore, 

expect  that  a  truncation  of  the'  matrices  B  and  is  permissible 

even  at  fairly  small  values  of  u.  However,  it  should  be 

understood,  that  because  of  the  peaks  which  occur  in  the  fund  ion 

<(u,r.)  for  u  close  to  1  ,  it  is  mandatory  that  the  integrations 

needed  in  the  determination  of  the  matrix  elements  of  B_  and  B . 

I  •> 

be  carried  out  with  sufficient  precision. 
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To  determine  quantitatively  to  what  extent  these  observations 
.ire  correct  ,  one  forms  the  eiqonvnlues  and  eigenvectors  of  the 
■  runcjtcd  Matrix,  B9  relative  to  H ^ .  The  elements  of  the  matrices 
a  r.  •  complex  because  the  functions  . ^  are  complex.  For  a  comparison 
■ '  f  the  e i qcnvcctors ,  one  uses  the  Hermit ian  scalar  nroduct .  If 


v  =  V  — 

1  a.  ’2  b. 

i  i 

t  hen  , 

f V x  ,  V2)  =  a.b.*  (89) 

v.licfc  the  asterisk  denotes  the  conjugate  complex.  Then  we  set 


cos 


!  tvr  V2]  I 
fVj,  V,]172  [V2 , 


(90) 


I  -  car:  ;  e  shown  that  1  =  0  only  if  =  const  V2 .  For  the  relative 
i ::  tin.  eigenvalues,  we  form 


k 


(91) 


A 
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SECTION  VIII 
NUMERICAL  RESULTS 


In  this  section  results  which  are  of  interest  for  practical 
work  are  collected. 

Figure  1  shows  a  number  of  curves  u  =  const.  They  can  be 
considered  as  ideal  boundaries, for  they  take  the  fact  into 
account  that  the  profile  (which  is  the  inner  boundary)  lies  at 
or  close  to  the  x  axis.  After  the  transformation  of  the  flow 
field  into  a  strip  of  the  r,n  plane  one  obtains  a  particularly 
simple  representation  of  the  potential  which  holds  throughout  the 
flow  field  and  displays  clearly  the  dominant  effects  at  large 
values  of  u  (or  O.  This  holds  for  the  underlying  steady  flow 
field,  as  well  as  for  superimposed  harmonic  perturbations. 
Deviations  from  these  boundaries  which  may  be  desirable  because 
of  practical  considerations  are  admissible.  In  fact,  the 
far  field  conditions  studied  here  hold  for  rather  arbitrary 
boundaries.  Which  of  these  boundaries  one  chooses  depends  upon 
the  character  of  the  partial  differential  equation.  In  deriving 
the  far  field  conditions  the  assumption  is  made,  that  in  the 
field  outside  of  this  boundary  one  deals,  after  the  necessary 
transformations  have  been  carried  out,  with  the  Helmholtz  equation. 

Part  of  the  transformations  which  lead  to  the  Helmholtz 
equation  is  a  Prandtl -Glauer t  transformation  in  the  physical 
plane.  The  ellipses  in  the  physical  plane  which  by  this  trans¬ 
formation  are  mapped  into  the  curve  u  =  5  are  shown  for  different 
Mach  numbers  in  Figure  2. 

The  constant  ;>  which  appears  in  the  Helmholtz  equation 
depends  upon  the  reduced  frequency  \\ '  -  L/u .  The  curves  of 
Figure  3  show  ..  versus  the  Mach  number  for  different  values  of 
a ' .  This  graph  allows  one  to  determine  the  range  of  values 
of  ..  which  is  of  practical  importance. 

It  was  m*--r.t  one!  that  an  iterative  solution  of  the  Helmholtz 
equation  is  possibb'  (but  perhaps  not  practical)  up  to  a  frequency 
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tor  which  a  solution  with  the  homogeneous  boundary  condition 
1  -  0  at  the  plate  4>  =  0  at  the  far  boundary  exists.  Figure  4 

shows  the  pertinent  values  of  1/ji  versus  u.  (The  frequency 
limitation  becomes  more  stringent  as  the  value  of  u  for  the  outer 
boundary  is  increased.) 

The  remaining  data  refer  to  the  far  field  conditions  directly. 
Let  us  first  repeat  the  leading  ideas.  The  far  field  conditions 
are  represented  by  an  infinite  system  of  linear  equations.  Two 
different  terms  of  this  system  have  been  considered,  namely  the 
representation  of  the  far  field  conditions  in  terms  of  products 
of  ordinary  and  radial  Mathieu  functions  in  the  p ,  Z,  system,  and 
the  representation  in  terms  of  products  of  Hankel  functions  of  r 
and  trigonometric  functions  of  6.  These  two  representations  are 
equivalent  to  each  other.  In  numerical  work  only  a  finite  number 
of  equations  of  these  systems  can  be  used.  The  equations  retained 
are,  of  course,  individually  correct;  the  approximation  lies  in 
the  fact  that  only  a  limited  number  of  these  conditions  are  satis¬ 
fied.  The  tar  field  conditions  do  not  introduce  errors,  but  they 
will  admit  certain  errors. 

The  partial  differential  equation  can  be  reduced  to  an  in¬ 
finite  system  of  ordinary  differential  equations.  Different  forms 
of  the  system  may  arise  if  different  representations  are  chosen. 

We  have  studied  a  representation  in  terms  of  Mathieu  functions  of 
and  radial  Mathieu  functions  of  C  alternatively  m  terms  of  a 
development  in  trigonometric  functions  with  respect  to  r,  with 
coefficients  that  depend  upon  u.  Here  errors  arise  by  the  trunca¬ 
tion  of  the  system.  We  have  studied  how  the  truncation  of  tfie 
system  representing  the  partial  differential  equation  interacts 
with  the  truncation  of  the  system  which  gives  the  far  field 
conditions. 

The  representation  of  the  flow  field  in  terms  of  Mathieu 
functions  converges  for  the  Helmholtz  equation  throughout  the 
field.  The  only  effect  of  the  truncation  of  the  system  is  seen 
in  the  boundary  conditions  at  the  plate  which  in  general  are 
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satisfied  only  approximately.  The  representation  in  terms  of 
Mathieu  function  clearly  expresses  the  directional  character 
of  the  radiation  (provided  that  the  inner  boundary  is  given  by  a 
slit  along  the  x  axis) .  Jn  account  of  these  observations  thi s  :<•- 
presentation  of  the  flow  field  can  be  regarded  as  best  from  a 
theoretical  point  of  view.  If  one  uses  this  representation  m 
a  practical  computation  then  the  flow  field  at  a  distance  can  be 
approximated  with  a  minimum  number  of  terms.  If  one  applies  to 
such  a  representation  a  formulation  of  the  far  field  conditions 
in  terms  of  Mathieu  function,  then  the  far  field  conditions  arc 
perfectly  satisfied. 

For  a  procedure  of  this  kind  one  needs,  of  course,  a  repre¬ 
sentation  of  the  Mathieu  functions.  Table  1  gives  for  =  0 .  a , 
y  =  2.0  and  y  =  3.5,  the  Fourier  coefficient  of  the  development 
of  (n)  defined  in  Fa.  (34),  for  the  five  lowest  eigenvalues. 

Let  us  consider  the  coefficients  a.  ,  as  elements  of  a  matrix, 

'  ;  K 

then  one  has  for  y  =  0  only  terms  in  the  main  diagonal  (foi  one 
deals  with  the  Laplace  equation).  As  y  increases,  the  important 
terms  are  found  in  the  vicinity  of  the  main  diagonal.  The 
"spreading  out"  of  the  Fourier  coefficients  does  not  increase  as 
one  proceeds  to  Mathieu  functions  with  larger  eigenvalues.  The 
presence  of  off-diagonal  terms  shows  the  directional  effect  of 
the  radiation.  It  is  present  even  for  the  lowest  eigenvalue, 
but  the  effect  is  pronounced  only  for  larger  values  of  ...  (For 
the  values  of  y  considered,  it  is  fairly  weak.)  The  eigenvalues 
are  shown  in  Tabic  1  underneath  in  an  extra  row.  They  must  be 

f 

used  to  compute  the  values  of  f^(u)/f^(u). 

The  counterpart  to  a  representation  of  the  flow  field  and  of 
the  far  field  conditions  in  terms  of  Mathieu  functions  is  a  for¬ 
mulation  in  the  physical  plane  in  terms  of  Uankel  functions  of 
the  radius  r  and  trigonometric  functions  of  the  angle  .  Such  a 
representation  would  be  practical,  if  the  inner  boundary  is  given 
by  a  circle.  If  the  inner  boundary  is  a  plate,  it  introduces  an 
undesirable  singularity  at  the  origin.  For  a  combination  of  a 
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a  representation  of  the  tlow  field  in  terms  of  triqonometric 
functions  oi  0  (this  is  the  important  part)  and  ot  the  far  field 
conditions  in  terms  of  llankei  triqonometric  functions,  the  far 
t  i-eLd  conditions  would  be  perfectly  satisfied. 

We  return  to  the  formulation  of  the  problem  in  the  u,  ' 
plane  and  consider  the  representation  of  the  flow  field  in  terms 
of  a  development  of  ‘  in  terms  of  a  truncated  Fourier  series  in 
with  coefficients  that  depend  upon  u.  Such  a  representation 
can  express  a  solution  in  terms  of  Mathieu  function  which  we 
consider  as  optimal  only  approximately.  This  fact  is  reflected 
in  the  expressions  for  the  flow  field  that  are  compatible  with  the 
far  field  conditions.  In  the  above  analysis  we  have  adopted  the 
following  characterization  of  the  far  field  conditions.  The 
functions  and  :  are  represented  as  linear  combinations  of  the 
relative  eigenfunctions  of  the  matrices  and  B^, ,  the  ratio  of 
the  coefficients  of  corresponding  terms  in  the  linear  combina¬ 
tions  tor  :  and  :  is  then  given  by  the  relative  eigenvalues. 

For  the  nontruncated  system  (in  any  representation)  the  eigen¬ 
functions  are  Mathieu  functions  of  q  the  relative  eigenvalues 
are  given  by  the  ratio  f(u)/f'(u)  obtained  from  the  radial  Mathieu 
functions.  Table  2  refers  to  a  representation  of  the  flow  field 
in  terms  of  trigonometric  functions  in  r  while  the  far  field  con¬ 
ditions  are  expressed  by  means  of  Mathieu  functions.  (We  have 
restricted  ourselves  to  solutions  which  are  symmetric  with  respect 
to  the  x  and  the  y  axis  and  then  restricted  ourselves  to  the  first 
five  terms,  the  last  term  represents  eight  cosine  waves  over  the 
whole  contour,  u  =  const.) .  While  the  Fourier  expansion  of  the 
Mathieu  function  is  obtained  by  the  eigenvectors  or  an  infinite 
matrix  one  now  deals  with  the  eigenvectors  of  the  same  matrix 
truncated  to  the  size  given  by  the  number  of  terms  in  the  flow 
field.  The  error  in  the  functions  admitted  arises  because  the 
eigenvectors  and  eigenvalues  are  not  the  same  in  the  two  cases. 


Table  2  gives  for  the  )j  =  0.5,  |i  =  2.0  and  =  3.5,  the 
angle  between  the  exact  and  the  approximate  eigenvectors  and  tin 
discrepancy  in  the  eigenvalues.  One  cannot  expect  perfect,  agree¬ 
ment  for  later  eigenfunctions,  because  the  truncation  of  the 
Fourier  series  suppresses  Fourier  terms  which  are  important  in 
the  exact  eigenfunctions.  The  table  shows  that  the  first  few 
eigenvectors  and  eigenvalues  are  indeed  very  close  to  the  exact 
values.  The  errors  increase  with  p.  The  eigenvalues  found  hero 
are  needed  in  order  to  determine  from  Eq .  (37)  the  ratio 

f^/ (u) /fR ' (u)  for  the  value  of  u  chosen  for  the  outer  bounuury. 

It  would  be  equally  justified  if  one  computes  the  values  of 
f^(u)/f^'(u)  for  the  infinite  matrices.  Then  no  error  in  the 
eigenvalues  would  be  encountered . 


If  one  uses  the  same  representation  for  the  flow  field  inn 
expresses  the  far  field  conditions  by  a  truncated  system  of  Hank*..; 
functions  in  r  and  trigonometric  functions  in  one-  will  obtain 
slightly  different  results,  because  of  the  difference  in  the  fur 
field  conditions.  Table  3  shows  for  different  values  of  .  and 
of  u  the  angle  of  corresponding  eigenvectors  for  the  truncated 
and  nontruncated  matrices  and  B ^  and  the  error  m  the  eigen¬ 
values.  In  this  case  the  evaluation  of  the  matrix  elements  is 
critical .  The  matrix  elements  are  integrals  over  periodic  analyt i 
functions  of  n.  In  this  case  the  trapezoidal  rule  ct  ives  excel  leni 
results,  but  the  integrands  are  rather  peaked  for  small  values 
of  u  and  therefore  the  integration  interval  must  be  chosen  sr.-.a 1 1 
enough.  Tables  3  have  been  computed  for  an  interval  of  integra¬ 
tion  tt/32  but  actually  for  the  values  of  u  considered  the  results 
are  the  same  if  one  uses  the  interval  .i/16 .  Serious  d  i  scrot-jr.rn-; 
arise  however  fo"  an  interval  ' /8 .  The  shortest,  wave  ler.  :  t  h 
in  the  representation  of  •;  by  a  Fourier  series  in  •  trace  t  ‘  -  ■ 

5  terms  is  '"/4.  An  interval  ~/8  would  give  just  2  joint:-  :  •  ,s 
full  wave.  For  small  values  of  u  the  peaked  ness  eg  t  V,.  :  •. 

is  increased  because  of  the  character  of  the  lianke  1  *  ,  : 

then  smaller  integration  intervals  will  probably  b<  m  e.  -sn. .  ry  . 
However,  u  =  4  is  probably  small  enough  for  all  practica' 
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purposes  (see  Fig.  1).  The  deviations  from  the  ideal  case  are 
significant  only  for  the  highest  eigenvalue. 

Table  4  shows  for  u  =  3.5  the  angle  between  the  eigenvec¬ 
tors  for  a  flow  field  represented  by  means  of  trigonometric  func¬ 
tions  in  u  with  the  far  field  conditions  expressed  either  by 
Mathieu  functions  (ordinary  and  radial)  or  by  Hankel  functions 
in  r  and  trigonometric  functions  in  0.  The  angles  are  extremely 
small;  for  smaller  values  of  y,  they  are  even  smaller) .  This 
shows  that  the  deviation  of  the  eigenvectors  from  their  ideal 
values  shown  in  Tables  2  and  3  is  primarily  due  to  the  discreti¬ 
zation  of  the  partial  differential  equation  and  not  to  the  form 
of  the  far  field  conditions. 

The  approximate  far  field  conditions  of  Bayliss  Gunzburger 
and  Turkel  are  mainly  of  interest  if  one  represents  the  flow  field 
in  terms  of  a  truncated  Fourier  development.  We  found  that  the 
e i qenf unc tu on  decomposition  which  one  finds  for  these  boundary 
conditions  is  the  same  is  if  one  imposes  the  far  field  conditions 
in  the  form  of  a  truncated  system  of  Mathieu  functions  and  these 
errors  are  already  shown  in  Table  2.  The  main  error  occurs  be¬ 
cause  of  the  falsification  of  the  ratio  f ^ (u) /f ^  1  (u)  .  Table  5 
shows  this  error  for  different  eigenfunctions,  different  values  of 
and  different  values  of  u.  As  expected  the  error  is  small  for 
sufficiently  largo  u  and  ...  Assume  for  instance  that  errors  in 
the  fourth  and  fifth  eigenvectors  are  unessential,  because  their 
contribution  to  the  solution  is  small  and  that  a  falsification 
oi  f/f'  for  the  third  eigenvector  by  .05  is  admissible,  then 
boundary  conditions  cannot  be  used  for  y  =  .5  and  u  as  large  as 
1  ■,  they  can  bo  applied  for  r  =  2  down  to  u  =  5  and  for  y  =  3.5 
<iown  to  u  -  3 . 
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SECTION  IX 

GENERAL  OBSERVATIONS 


One  will  ask  to  what  extent  the  results  of  the  present 
analysis  can  be  carried  over  to  more  realistic  problems. 

The  transformations  which  in  the  present  case  lead  to  the 
Helmholtz  equation  can  be  useful  for  higher  values  of  ;j ,  because 
it  reduces  the  waviness  due  to  waves  that  travel  upstream. 

The  choice  of  the  boundary  of  the  computed  flow  field  as  a 
line  u  =  const  arises  because  of  the  shape  of  the  inner  boundary 
(here  the  plate) .  It  can  be  expected  to  be  advantageous  for  all 
two-dimensional  airfoil  problems,  but  some  modifications  are 
admissible . 

The  far  field  (especially  its  directional  characteristics) 
is  best  represented  in  terms  of  Mathieu  functions.  This  is  a 
consequence  of  the  shape  of  the  inner  boundary. 

The  idea  of  representing  <P  in  terms  of  Mathieu  functions 
throughout  the  flow  field  requires,  of  course,  that  one  employs 
coordinates  (or  a  mesh)  given  by  the  u,  n- system.  The  idea  of 
conformal  mapping  has  been  used  even  in  the  transonic  case,  there¬ 
fore,  a  precedent  to  this  procedure  exists.  The  use  of  a  Mathieu 
function  development  (as  well  as  Fourier  series  development  in  r) 
may  offer  difficulties  in  the  vicinity  of  the  air  foil,  where  the 
function  <j>  may  be  rather  peaked  (because  the  character  of  the 
boundary  conditions,  and  also  because  of  the  discontinuities  which 
may  be  encountered  in  a  transonic  flow  field) . 

The  representatives  of  far  field  conditions  in  terms  of 
Mathieu  functions  (even  if  the  boundary  is  not  exactly  a  line 
u  =  const) is  quite  feasible  (it  gives  best  accuracy  with  a 
minimum  number  of  terms) . 

The  representation  of  the  flow  field  in  terms  of  a  Fourier 
development  in  n  requires  a  few  more  terms  than  a  representation 
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in  terms  of  Mathieu  functions,  if  one  wants  to  capture  the  direc¬ 
tional  characteristics  of  the  far  field  equally  well.  We  men¬ 
tioned  above  that  in  the  representations  of  the  field  close  to 
the  profile  both  forms  have  the  same  limitations.  One  expects 
ti.at  other  representations  of  the  flow  field  will  have  the  same 
requirements  for  the  representation  of  the  distant  field  as  a 
Fourier  representation.  They  may  have  advantages  in  the  vicinity 
of  the  profile. 

Errors  in  satisfying  at  the  far  field  conditions  originate 
mainly  by  the  truncation  needed  in  the  partial  differential 
equation.  Whether  one  represents  the  far  field  conditions  by 
means  of  Mathieu  functions  or  by  Hankel  functions  in  r  and  tri¬ 
gonometric  functions  0  has  practically  no  effect,  if  the  outer 
edge  of  the  computed  parts  of  the  profile  is  given  by  a  line 
u  :•  4. 

The  approximate  far  field  conditions  of  Bayliss,  Gunzburger, 
and  Turkel  can  be  applied  only  if  y  and  u  are  sufficiently  large. 
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APPENDIX 
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LEGENDA  TO  FIGURES 


Figure  i 
Figure  2 

Figure  3 
Figure  4 


Curves  u  =  const  in  the  x,y  plane  after  a  Prandtl- 
Glauert  Coordinate  Distortion  has  been  carried  out. 

Lines  u  =  5  in  the  x,y  plane  for  different  Mach 
numbers,  if  the  Prandtl-Glauert  Coordinate  Distortion 
is  not  carried  out.  Such  curves  constitute  ideal 
outer  boundaries  for  a  computed  flow  field. 


Values  of  u  for  different  reduced  frequencies  u 
function  of  the  Mach  number. 


as 


Values  of  u  for  which,  at  a  given  value  (1/u)  an 
iterative  procedure  will  theoretically  converge  (the 
admissible  region  has  been  shaded) .  It  is  assumed 
that  in  the  iterative  procedure  the  values  of  4  at 
the  outer  boundary  will  be  recomputed  in  each  iteration 
step  and  then  kept  fixed  during  the  flow  field 
computation.  The  theoretical  limit  arises  for  a 
chosen  value  of  u  at  the  lowest  frequency  u  for  which 
a  standing  wave  flow  with  <1  =  0  at  the  outer  boundary 
and  r-  =  0  at  the  inner  boundary  can  arise. 


-«.M  -S.M  -Ml  -3  -Z.M  -J.M  l.»  ».••  B.**  3.H  4.m  S.H  «.#• 

Figure  1.  Constant  u  Contours  in  ,  y  Plane  With 


Pranctl-Glaueit  Coordinate  Distortion. 


F inure  2.  riach  'Cumber  Contours  for  u  =  5  in  x,  y 
Without  Prandt  1  -Clauer t  Coordinate  Distortion. 
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LEGKMDA  TO  TABLES 


a!  - 1 1  ■  1 


Table  2 


Table  3 


1 


[ 

i 


Tabic  4 


Tabic  5 


Fourier  I’oef  f ic i on t  s  of  *  i  vc  Mithieu  i'unet  ion::  wi :  i  !. 
a  re  symmetric  to  the  >:  and  y  Axes  and  e<  rr-.-sj-wn-i  i  nu 
iaqenvalucs. 

A  characteri zation  of  the  Far  Field  errors  for  a 
Flow  Field  expressed  by  Trigonometric  Functions  in 
■  and  Far  Field  Conditions  formulated  by  means  of 
Mathieu  Functions. 

The  table  shows  the  angle  between  the  relative 
eigenvector  of  the  two  matrices  which  appear  in  the 
far  field  conditions  computed  once  for  a  perfect 
representation  of  the  flow  field  and  a  second  time 
for  a  Fourier  scries  truncated  to  five  terms.  In 
practice  these  eigenvectors  are  computed  from  the 
matrix  which  give  a  Fourier  representation  for 
Mathieu  functions.  The  infinite  matrix  has  been  approx- 
mated  by  a  nine  by  nine  matrix,  the  approximate  eigen¬ 
vectors  arise  from  a  matrix  truncated  to  five  by  five. 

Also  given  are  the  relative  errors  in  the  eigenvalues 
for  these  matrices.  These  eigenvalues  serve  to  deter¬ 
mine  the  value  of  f(u)/f'(u)  from  F.q .  (35).  The  values 

f(u)/f'(u)  are  the  eigenvalues  which  appear  in  the  charac¬ 
terization  of  the  far  field  conditions. 


A  Character i zation  of  the  Far  Field  Errors  for  a  Flow 
Field  expressed  by  Trigonometric  Functions  in  ■  and 
Far  Field  Conditions  given  by  Henkel  Functions  in  the 
Radius  and  Trigonometric  Functions  in  the  Angle  in  the 
Physical  Plane. 

The  table  shows  the  angle  between  the  relative  eigen¬ 
vectors  of  the  two  matrices  which  appear  in  the  far 
field  conditions  computed  once  for  the  perfect 
representation  of  the  flow  and  a  second  time  if  the 
Fourier  series  is  truncated  to  five  terms.  The  error 
in  the  eigenvalues  refers  to  the  matrix  which  appears 
in  the  far  field  conditions. 

Angle  between  the  Relative  Eigenvalues  between  the 
Two  Matrices  which  appear  in  the  Far  Field  Conditions 
for  a  Flow  Field  represented  by  a  Truncated  Series  in 
Trigonometric  Functions  in  r  .  The  far  field  conditions 
are  expressed  once  in  terms  of  Mathieu  functions  and 
a  second  time  in  terms  of  Hankel  functions  in  the 
radius  and  trigonometric  functions  in  the  angle. 

Relative  Error  in  the  Eigenvalues  for  the  conditions 
of  Bayliss,  Gungburger  and  Turkel ,  for  a  Flow  Field 
Represented  by  Trigonometric  Functions  in  The  Error 

in  the  Eigenfunctions  is  the  Same  as  in  Figure  2. 
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*See  Legenda  on  page  59.  j 


r 


TAB !  ,E  2* 
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See  Leqenda  on  page  59. 
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